Strong eigenfunction correlations near the Anderson localization transition. 
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We study overlap of two different eigenfunctions as com- 
pared with self-overlap in the framework of an infinite- 
dimensional version of the disordered tight-binding model. 
Despite a very sparse structure of the eigenstates in the vicin- 
ity of Anderson transition their mutual overlap is still found 
to be of the same order as self-overlap as long as energy sepa- 
ration is smaller than a critical value. The latter fact explains 
robustness of the Wigner-Dyson level statistics everywhere in 
the phase of extended states. The same picture is expected 
to hold for usual d— dimensional conductors, ensuring the s^ 
form of the level repulsion at critical point. 

Recently, there was considerable growth of interest 
towards constructing a unified picture of wave func- 
tions and energy levels statistics for disordered conduc- 
tors in the vicinity of Anderson metal-insulator transi- 
tion |p]-[T^. Coming from the metallic phase, a typical 
wavefunction ^'^(r) is extended and covers all the sam- 
ple volume randomly, but uniformly. When system ap- 
proaches the point of Anderson transition Ec, these ex- 
tended eigenfunctions become less and less homogeneous 
in space showing regions with larger and smaller ampli- 
tudes and eventually forming a multifractal structure in 
the vicinity of Ec- 

To characterize the degree of non-homogeneity quan- 
titatively it is convenient to use the inverse participation 
ratio (IPR) I{E) = J dra{r,E), where 



a(r, E) = {\^dr)\')E ^ A ( ^ |vl'„(r)|4<5(£ - E, 



(1) 



A is the mean level spacing and the angular brackets 
stand for the disorder averaging. For extended states this 
quantity is inversely proportional to the system volume: 
I{E) = C{E)L^'^, with L and d standing for the system 
size and spatial dimension, respectively. The coefficient 
C in this relation measures a fraction of the system vol- 
ume where eigenfunction is appreciably non-zero. For 
random homogeneous states C ~ 1, whereas close to the 
mobility edge E — Ec it becomes large and diverges like 
C{E) (X \E — Ec\~^, ^ > signalizing of increas- 
ing sparsity of eigenfunctions. Just at the mobility edge 
eigenfunctions occupy a vanishing fraction of the system 
volume and IPR scales like I{E) cx L~'^^^, r] > 0. Such 
a behavior reflects fractal (actually, multifractal p^,^,p[ ) 
structure of critical eigenstates. At last, in the insulating 
phase any eigenstate is concentrated in a domain of finite 
extension ^; and IPR stays finite in the limit of infinite 



system size L —>■ oo. 

This transparent picture serves as a basis for qualita- 
tive understanding of spectral properties of disordered 
conductors. Indeed, as long as eigenstates are well ex- 
tended they overlap substantially and corresponding en- 
ergy levels repel each other in the same way as do eigen- 
values of large random matrices studied by Wigner and 
Dyson. As a result, the Wigner-Dyson (WD) statistics 
describes well energy levels in a good metal jl^Jl^. In 
contrast, in the insulating phase different eigenfunctions 
corresponding to levels close in energy are localized far 
apart from one another and their overlap is negligible. 
This is the reason for absence of correlations of energy 
levels in this regime - the so-called Poisson statistics. 

However, close to the transition point such a reasoning 
should be used with caution. Naively one may expect 
that sparse (multifractal in the critical point) eigenstates 
fail to overlap, that would result in essential weakening 
of level correlations close to the mobility edge and van- 
ishing level repulsion aX E — E,.. However, a thorough 
investigation shows 0:^,^^^ that even at the mobility 
edge levels repel each other strongly, though the whole 
statistics is different from the WD one. One of the main 
purposes of the present paper is to resolve this apparent 
contradiction. We will show that critical eigenstates for 
nearby levels are strongly correlated and overlap well in 
spite of their sparse structure. 

The overlap of two different eigenstates VPi and 5" j cor- 
responding to energy separation \Ei — Ej\ = uj can be 
characterized by comparing the correlation function 



a{r,E,oj) = {\^,{v)\^\^,{v)\^) 



\^,ir)\^\^,ir)\HiE E,)6{E + - E,) ) (2) 



at UJ ^ with a(r, E), the latter function playing in such 
a context the role of the eigenfunction self-overlap. Here 
R2{uj) denotes the two- level correlation function. 



R2{uj) = - E{)5{E + uj~Ej) 



(3) 



To study the function a{r, E, uj) analytically, we consider 
an exactly solvable model of the Anderson transition - so- 
called sparse random matrix (SRM) model Jl7| , p^ . This 
model deals with large N x N matrices (real symmetric 
or Hermitian) whose entries Hij are independent random 
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variables characterized by the following probability dis- 
tribution: 



V[H,,) = {l-p/N)5{H,,) + (p/N) h{H,,) 



(4) 



where h{z) = h{—z) is any even distribution function 
having finite second moment. The parameter p > has 
a meaning of mean nonzero elements per row. For p > 1 
the model describes a connected graph having locally a 
tree-like structure with random connectivity. On a larger 
scale there exist large loops involving typically of the or- 
der of In N sites. In the limit — > oo the influence of 
loops is negligible and the model belongs to the same 
universality class as disordered tight-binding model on 
the infinite tree (Bethe lattice). The latter model was 
thoroughly investigated by various methods [|l9|,0 and 
shown to possess the Anderson transition. 

Additional interest to structure of eigenfunctions on 
tree-like structures is attributed by a recent paper 
There, the problem of quasi-particle lifetime induced by 
Coulomb interaction in mesoscopic samples was mapped 
onto a disordered tree-like tight binding model in the 
Fock space, which is similar to the Bethe lattice and SRM 
models and thus undergoes the localization transition. 
The fact that in the vicinity of the transition the cor- 
responding extended eigenfunctions are extremely sparse 
made the authors to conclude that the level statistics dif- 
fers from the WD one for such a regime because "eigen- 
functions do not talk to each other" ||2l[ . Such a conclu- 
sion would be, however, at variance with the results of 
explicit calculation of level-level correlation function per- 
formed in the framework of SRM model , where it was 
shown that i?2(w) is given by the WD form everywhere 
in the region of delocalized states up to the transition 
point E — Be- 
lt is necessary to mention that SRM model has some 
considerable advantage when compared with more con- 
ventional Bethe lattice model. Namely, quantities like 
IPR are not unambiguously defined on the Bethe lattice 
in the phase of extended states. Indeed, they require the 
consideration of a large but finite lattice to be well de- 
fined and their limiting behavior may depend crucially on 
boundary conditions imposed In contrast, all sites 
of SRM model are essentially equivalent and the model is 
free from boundary problems. The general expression for 
IPR was derived in 1 1^ and its critical behavior analyzed 
in(2|]. 

Actually, the SRM model can be used to construct an 
effective mean-field theory of Anderson localization ||2^ ] 
valid at d = oo [Q. Critical properties of such theo- 
ries first discovered on the level of non-linear ct— models 
p5| , p6t turn out to be quite peculiar. In particular, the 
coefficient C{E) diverges close to the transition point 
like C{\E - Ec\ < E^) oc exp (const|£; - Ec\-^/^) || in 
contrast to expected power-law behavior for conventional 
d— dimensional systems. The origin of such a critical de- 



pendence was explained in 24 1 and stems from the 
fact that C (E) is determined essentially by the " correla- 
tion volume" V{£,) (i.e. number of sites at a distance 



smaller than correlation length which is exponen- 
tially large, V{£,) oc exp(const^), for tree-like structures, 
whereas V{^) oc ^'^ for a c?— dimensional lattice. Having 
this difference in mind, one can translate all the results 
obtained in the framework oi d — oo models to their 
finite-dimensional counterparts ||2^ . 

To calculate the overlap function defined in Eq.(|^) 
we follow |27) and use the identity relating a{r, E) and 
a{r^E,u}) to advanced and retarded Green functions 



E±iri-Ei ' 



27r' 



A-^a{r, E)6{uj) + A'^ R2{Lu)a{r, E, uj) 



(5) 



= Re [(G^(r, E)G^{r, E + uj) - G^{r, E)G'^\r, E + uj))] 

where i?2(w) is non-singular part of the level- level cor- 
relation function: R2{'jj) — i?2('^) + (5(w/A). Let us 
consider for definiteness the ensemble of real symmet- 
ric SRM, corresponding to systems with unbroken time- 
reversal invariance. For any site index r = l,...,N 
we introduce one eight-component supervector <I>t = 

consisting of two four-component supervectors 

= [(i>aM,4>a,b2-,(t>oj-, -<l>aj^-, whcrc iudiccs a = R,A 
and h, f are used to label advanced-retarded and boson- 
fermion subspaces, respectively. The ensemble-averaged 
products {C^G"') for RSM model in the Hmit iV > 1 
can be extracted from the paper and the Appendix 
D of the paper E3] and is given by: 



G"{r,E)G" ir,E + Lu)j= i^l 

DQ {4>a,bl4><y,bl4>c!' ,bl4'a' ,bl) 

= / . .) exp 



exp ( ^^^^^ Sir QA 



(6) 



The function 5t($) = 5o(*^rtr$; $ti$) satisfies the 
integral equation: 



5t(*) = {[hp ($tL*) -1]) 



St 



(7) 



where hF{t) — j dze^^'^^h{z) is the Fourier-transform of 
the distribution of nonzero elements of the SRM. The 
8x8 supermatrices T satisfy the condition T''LT — L 
where L = diag(l, 1, 1, 1, — I, — 1, 1, 1) and belong to a 
graded coset space whose explicit parametrisation can 
be found in ]l5| , p8t . The supermatrices Q are expressed 
in terms of T as Q = T^^LT. At last, the matrix A — 
diag(l, 1, 1, 1, —1, —1, —1, —1), and the density of states 
p is expressed in terms of the solution of the equation 
Eq.(0) as p{E) = -2go^/(7rB2), where B2 = J dzh{z)z^ 
and go,, = dgo{x, y)/dx\^,y^o] x = y = 

When deriving Eq.(H), evaluation of a functional inte- 
gral by the saddle-point method has been employed, see 
An accurate consideration shows that 
ure is legitimate as long as: i) the matrix 



details in 
such a procec 
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size N (playing in our model the role of the volume) is 
large enough (much larger than the coefficient C{E) de- 
termining the size dependence of IPR, see above); and 
ii) the energy difference lo is small enough (much smaller 
than C^^{E)). Though C{E) is exponentially large near 
the transition point, it depends on the energy E only, 
so that when we keep E fixed and increase the system 
size N, the number of levels in the interval C~^{E) gets 
arbitrarily large, since the level spacing scales as 1/N. 

Expanding both sides of Eq. ^ over "if one can express 
{4'a,bi^'y,bi4>a' ,bi4>a' ,bi) hi tcrms of the matrix Q as 



'cr,61' 

4! 



■,bl'. 
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'90-xx h\ ~^ 9Q,xyQbibi ~^ 90,'t 



{4'RM4>R,bl4>A,bl't'A,bl) gj, 

= ^"5- [90.XX {QbibiQbibi + '^QbibiQ 

90,xy + 90,yy\ : 



AR \ , 



(8) 



where 50, = d'^90/ dx^\x,y=o] go,, 
9o.,xy = d'^9o/dxdy\x,y=ii, and B4 



= d^go/dy^\x,y=o; 
J dzh{z)z^. This 



allows us to represent right-hand side of Eq.(H) in the 
following form: 



2tt^ A~ia(r, E)5{uj) -I- /^'^ R{ijj)(3{r, E, uj) 
= --^9oxx^e {{Q^f Q^^ + 2Q^_^ Q^^ 
1 \nRR nRR ^ n^A n^A - 

9 \yibibi^bibi + ^bibi^bibi. 



(9) 



where 



(••■)q 



dg(...)exp(^!^Str QA 



The integrals over Q— matrices are the standard ones [ p^ 
yielding: 



(QbibiQbibi) o - ^ 



2iA 



/nRR nRR \ — / r^AA 

\Wbibi^bibi/c 



7r(u;-hz0)' 
(Qbi^iQbiii) Q 
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(10) 



where r'^'^ {lu / A) is the level correlation function in the 
Gaussian Orthogonal Ensemble. Substituting this in 
Eq.(||), we finally find 



a{r,E,Lu) 



1 / 1 4g0,a:x 

—a(r,E) — — ^ ^ 



(11) 



The coefficient 1 /3 in Eq. (|ll|) corresponds to the case 
of unbroken time reversal symmetry (orthogonal ensem- 
ble) . For the unitary ensemble (broken time reversal sym- 
metry) the same consideration yields the coefficient 1/2 
instead, so that the general relation reads 



cr(r, E,u) 



(3 + 2 



a{r, E) 



(12) 



where (3 is the conventional symmetry parameter equal to 
/3 = 1 (2) for the orthogonal (resp. unitary) ensembles. 
This relation between the overlap of two different eigen- 
functions (7{r,E,uj) and self-overlap a{r,E) constitutes 
the main result of the present publication. It is valid 
everywhere in the phase of extended eigenstates, up to 
the mobility edge E = Ec, provided the number of sites 
(the system volume) exceeds the correlation volume. In 
particular, it is valid in the critical region \E — Ed <C E^ 
where a typical eigenfunction is very sparse and self- 
overlap (hence, IPR) grows like exp (const|_E — Ec\^^^^) 

ii- 

Eq.(p^ implies the following structure of cigcnfunc- 
tions within an energy interval SE = uj < C~^{E). 
Each eigenstate can be represented as a product Vf^ (r) = 
Eif) ■ The function $_E(r) is an eigenfunction en- 
velope of "bumps and dips" which is smooth on a micro- 
scopic scale comparable with lattice constant. It is the 
same for all eigenstates around energy E, reflects under- 
lying gross (multifractal) spatial structure and governs 
the divergence of self-overlap at critical point. In con- 
trast, "01 (^) is Gaussian white- noise component fluctuat- 
ing in space on the scale of lattice constant. It fills in 
the " smooth" component e {r) in an individual way for 
each eigenfunction, but is not critical, i.e. is not sensitive 
to the vicinity of the Anderson transition. These Gaus- 
sian fluctuations are responsible for the factor (3/ {(3 + 2) 
(which is the same as in the corresponding Gaussian En- 
semble) in Eq.(^. 

As was already mentioned, this picture is valid in 
the energy window SE C^^{E) around the energy 
E\ the number of levels in this window being large as 
SE/A ^ NC~^(E) 1 in the thermodynamic hmit 
N ^ 00. These states form a kind of Gaussian Ensem- 
ble on a spatially non-uniform (multifractal for E —^ Ec) 
background $£;(r). Since the eigenfunction correlations 
are described by the formula (|l^) , which has exactly the 
same form as in the Gaussian Ensemble, it is not surpris- 
ing that the level statistics has the WD form everywhere 
in the extended phase [0 . 

We believe on physical grounds that the same picture 
should hold for a conventional d-dimensional conductor. 
First of all, the general mechanism of the transition is 
the same in d < 00 and d — 00 models. Furthermore, the 
sparsity (multifractality) of eigenstates near the transi- 
tion point takes its extreme form for d = go models [ p4[ , 
so that since the strong correlations ( |l^ ) take place at 
d = 00 it would be very surprising if they do not hold at 
finite d as well. Finally, Eq. (|l^ ) was proven by an explicit 
calculation in the weak localization regime p^ , where 
a{r,E,oj) = ■^a{r,E) = V-^[l + U{r,r)], with V be- 
ing the system volume and n(r, r) the diffusion propa- 
gator. 

Replacing C{E) by the d-dimensional correlation vol- 
ume ~ we conclude that for E close to Ec eq.dl2) 
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should be valid for ui < A^, where ex 1/^'^ is 
the level spacing in the correlation volume. For larger 
uj, (7{r,E,uj) is expected to decrease as cj^''/"^' accord- 
ing to the scaling arguments pP,]!^, so that we find 
a{r,E,uj)/a{r,E) ~ (tj/A^)"''", up to a numerical co- 
efficient of order of unity. Again, for any value of the 
energy E in the delocalized phase, taking the system size 
L large enough, L 3> ^, we have a large number of levels 
SE/A ^ A^/A oc {L/^Y in the energy window 5E where 
eq.(|IJ) holds, so that the level correlation will be of the 
WD form. 

Finally, let us consider what happens when we go from 
the critical regime large, but L 3> ■C) to the critical 
point ^ L). For this purpose, let us keep the system 
size L fixed and change the energy toward Ec, so that ^ 
increases. When ^ gets comparable to the system size, 
^ ^ L, we have Aj ~ A. This is the border of appli- 
cability of the above consideration. Correspondingly, we 
find 
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(7{r,E,uj)/a{r,E) - 1 



< A 



(13) 



and a{r,E,uj)/a{r,E) - (w/A)"''/'* for u; > A. When 
E approaches further Ec, the correlation length ^ ^ i 
gets irrelevant, so that these results will hold in the crit- 
ical point (^ = oo). Of course. Eg. (p^) is not sufficient 
to ensure the WD statistics in the critical point, since 
there is only of order of one level within its validity range 
SE ~ A. Indeed, the numerical simulations show that 
the level statistics on the mobility edge is different from 
the WD one 0E|-§- 

However, eq.(|l^) allows us to make an important con- 
clusion concerning the behavior of i?2(w) at small lo < A, 
or, which is essentially the same, the behavior of the 
nearest neighbor spacing distribution P(s), s = uj / A, at 
s < 1. For this purpose, it is enough to consider only two 
neighboring levels. Let their energy difference be wq ~ A. 
Let us now perturb the system by a random potential 
V{r) with {V{r)) = 0, {V{r)V{r')) = T5{r-r'). For the 
two-level system it reduces to a 2 x 2 matrix {Vy}, i,j = 
1,2, with elements V^j = / d''rV^(r)**(r)^'j (r). The 
crucial point is that the variances of the diagonal and off- 
diagonal matrix elements are according to Eq.([l3|) equal 
to each other up to a factor of order of unity: 



{V,\)/{\V,l\) ^<j{r, E,^)/a{r,E) ^ 1 



(14) 



The distance between the perturbed levels is given by 
UJ = [{Vii - V22 + ujo)"^ + |142p]^/^- Choosing the am- 
plitude of the potential in such a way that the typical 
energy shift Vn ~ A and using eq.(|lj) and the stan- 
dard symmetry consideration, we immediately conclude 
that in the critical point P(s) ~ cps^ for s <^ 1 with 
a coefficient c^g of order of unity, in agreement with the 
numerical findings [||J|,01 . 
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